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ABSTRACT 

We  are  concerned  with  the  application  of  two  finite 
difference  schemes  of  second  order  accuracy  to  the  equations 
of  viscid  and  Invlscld  flow.   The  first  Is  the  Lax-Wendroff 
scheme,  the  second  an  Iterative  scheme.  A  stability  criterion 
Is  obtained  for  the  Iterative  scheme  by  the  method  of 
von  Neumann.   An  empirical  stability  criterion  Is  obtained 
for  the  Lax-Wendroff  scheme  as  applied  to  viscid  flow.   The 
accuracy  of  these  methods  and  their  effectiveness  for  flows 
which  contain  a  shock  are  also  discussed. 
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NYO-9188 

ON  CERTAIN  FINITE  DIFFERENCE  SCHEMES  FOR 
THE  EQUATIONS  OF  HYDRODYNAMICS 

1 .   Introduction . 

We  will  consider  two  finite  difference  schemes  for  the 
solution  of  the  Navier-Stokes  equations  in  one  space  di- 
mension.   The  first  is  a  method  used  by  Lax  and  Wendroff  [l], 
The  second  method  is  based  on  an  approximate  solution  of  the 
centered,  implicit  difference  equations.   Since  this  approxi- 
mate solution  is  obtained  by  successive  substitutions,  we 
will  refer  to  this  second  method  as  an  iterative  method. 

By  using  the  method  of  von  Neumann  we  will  obtain  a 
stability  criterion  for  the  iterative  method  as  applied  to 
the  equations  for  inviscid  flow.   This  method  is  stable  if 
3,  4,  7j  8,  ...    iterations  are  used  per  time-step  and  un- 
stable if   1,  2,  5,  6,     ...     iterations  are  used.   We  will 
obtain  an  empirical  stability  criterion  for  the  Lax -Wendroff 
method  as  applied  to  the  Navier-Stokes  equations . 

Both  of  these  methods  have  a  truncation  error  of  second 
order.   We  will  describe  the  results  of  calculations  in 
which  both  of  these  methods  were  applied  to  a  problem  for 
which  an  exact  solution  is  known.   Thus  we  will  obtain  an 
indication  of  the  accuracy  of  these  methods. 

P.  Lax  has  applied  the  Lax -Wendroff  method,  as  well  as 
a  method  of  first  order  accuracy,    to   inviscid  flows  which 
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contain  a  shock  of  moderate  strength  [l,2].   He  used  the 
equations  of  motion  In  conservation  form  In  order  to  handle 
the  discontinuity  without  the  Introduction  of  artificial 
viscosity  or  the  use  of  shock  fitting.   We  will  compare 
the  results  of  computing  with  the  Iteration  method  and  the 
Lax-Wendroff  method  for  Invlscld  flows  which  contain  a  shock. 
It  will  be  shown  that  the  Iteration  method  with  the  equations 
In  conservation  form  does  not  work  as  well  as  the  Lax-Wendroff 
method  where  the  equations  are  not  In  conservation  form. 

2 .   The  Finite  Difference  Equations . 

We  let   p,  p,  u,  e,  m,  j,      and  R  denote  density,  pressure, 
velocity,  energy,  momentum,  ratio  of  specific  heats,  and 
Reynolds  number   (m  =  pu  and   e  =  p-  pu   +  —^)  •      The  Navler- 
Stokes  equations  can  be  written  as 


where 


V  = 


m 


^  +  ^  +  M.  0 
^  t    ^x    ^x 


f(v)  = 


in 


7me  _  7-1   m 
P      2    p^ 


3 


(^_i)e  -  z^ni 


2 
2   P  J 


g  = 


m 


_  4_  m  d(p 
3R  p  ^x 


m 


_  K_   a(p) 

5R  ^x 


We  denote  the  forward  and  backward  differences  by 

V   =  v(t,x  +  zlx)  -  v(t,x)   and   v.  =  v(t,x)  -  v(t,x  -  Ax) 

X  X 

and  the  centered  difference  by  v.v  =  /2(v.-  +  v  ).   We  use 


X 


X' 


the  notation  v.  =  v(t  ,x.)   and  denote  the  Jacoblan  of 
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f(v)   with  respect  to   v  by  A.   The  Lax-Wendroff  method  is 
based  on  the  expansion 


(2.1) 


a 


V 


n 


at 


n 


iH 


+ 


a  g 


bxbt/ . 
/  1 


hf   ^  A   n 

(note  that  [•^—)      -    (A  -s— )  )  .   In  order  to  obtain  a  difference 


''^ 


^ 


approximation  for  the  term   ^  g/Sx^t  we  need  a  difference 
approximation  for  du/dt.   This  Is  obtained  from  the  Navler- 
Stokes  equation 


au 

at 


-  u 


au  _  _i  a^  ^    4     a  u 
ax     p  ax     3Rp    ax 


by  replacing  derivatives  by  difference  quotients.   If  we 
denote  this  approximation  by  A,u  and  denote  the  similar 
approximation  for  ag/at   by  A,g,   then 


A^u  =  -  uu^  -  1  p^  +  Jt 


3Rp 


XX    and 


(A^g)*  =  fo,  -  —  (u(A^u)^  +  (A^u)u^,  -  —  (A^u)^ 
N     3R  3R 


(g  denotes 
the  trans- 
pose of  g) 


The  Lax-Wendroff  method  Is  obtained  by  replacing  derivatives  with 
difference  quotients .   The  difference  equations  are 


n+1    n 

V     =  V 


A 


/„n  ,   n  \    A 


2  r     n        n 
(AAv^)_  +  (A^g)_ 

X  X 


where 


At 


A  =  ^  ,  Av,  =  -  f 
Ax  '    t       X 


4 


g^,   and  g  =  (0,  -  ^  uu^. 


3R  X 
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In  order  that  the  truncation  error  be   O(Ax^)  the 
differences   f/\   and  g/\  must  be  centered.   Since  the  last 
term  on  the  right  Is  of  second  order,  the  differences  In  this 

term  need  not  be  centered.   In  the  calculations  the  boundary 

n        n 
values  of   v,   that  Is   v-,   and  v^,   are  known  and  the  values 

n 
of  V.   for   1  =  2,3,...M-1   and  n  =  2.,^,.,.        are  computed. 

Since  the  second  order  term  In  the  difference  equations  requires 

,,     -1       nnnn       ^n,        ,     n+1 
the  values   v._p,  v.  ,,  v.,  v.  , ,   and  v.  „    to  compute   v.   , 

we  must  modify  this  term  when   1=2   or   1  =  M-1 .   The  modi- 
fication merely  requires  that  certain  forward  differences  be 
replaced  by  backward  differences,  or  conversely. 

We  also  performed  calculations  using  the  Lax-Wendroff  method 
for  the  Navler-Stokes  equations  where  the  latter  are  not  written 
In  conservation  form.   In  this  case  the  equations  are 

3  p      ^u    dp 

—^  -   -   p —  -  u— ^ 

dt      dx    Sx 


(2.2)   ^=  -  u^  -  7P—  +  ^^^^^^(^)^ 
dt      dx     dx     3R   dx 


du  ^  _  ^^  _  1  ^  ^.   ^   d^u 

dt     dx   p  dx   3RP  ^y? 

The  difference  equations  are  obtained  In  the  same  manner  as  before. 

The  Iteration  method  Is  based  on  an  Implicit  form  of  the 
difference  equations .   These  are 

n+1 
v  =  V      - 


n        A    /pn+1    .    ^:,n  \         A    An+l   ^  -n\ 
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where   g  (0, [-) ^7;^-)  _  . (-)  _ 

5RAx   ^  X   3RAx  ^      ^   XX  3RAx  ^   xx 


The  Iteration  method  is  defined  by 


o    n 
w   =  V 


n+1    p 
V    =  w^ 


w^+1  =  v'^  -  ^  ((f(w^))^  +  (f(v'^))^) 


-  ^  (i(w^^+  g(v^)A). 


This  merely  amounts  to  an  approximate  solution  of  the 
Implicit  difference  equations  by  successive  substitutions 
where   p   is  the  number  of  Iterations .   We  use  the  same  value 
of   p   at  all  time  steps,  i.e.  for  all  values  of  n.   The 
truncation  error  of  this  method  is   OCAx''^)   provided  p  ^  2. 
Since  the  implicit  difference  equations  are  unconditionally 
stable,  we  might  conjecture  that  this  iteration  method  has 
a  favorable  stability  criterion.   Unfortunately  this  is  not 
the  case.   We  will  show  that  the  method  is  unconditionally 
unstable  if   p  =  1,  2,    5,  6,  9,    10,     ...   and  conditionally 
stable  otherwise  with  the  condition  that   A   be  less  than 
twice  the  bound  that  Courant,  Priedrichs,  and  Lewy  discovered  [3] 

5 .    Description  of  the  Computed  Solutions. 

Professor  H.  Keller  suggested  the  use  of  this  problem. 
Ludloff  and  Filler  have  computed  the  solution  to  this  problem 
using  difference  schemes  of  first  order  accuracy  [6]. 
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We  compute  the  solution  between  x  =  0  and  x  -  c   at 

the  points   (t  ,x. )   where   x.  =  lAx,  t   =  nAt,   1=1,  2,...  M, 
f  n'  1  1      '   n      '       >      }  > 

n  =  1,  2,...    N,   MAX  =  c  >  1 .   We  assign  the  Initial  values  of 
u  as  follows   (0  <  d,  d  +  1  <  c) 

+  1 

u(0,x)  =   <^  0ix)    (arbitrary) 

0 

We  choose   0(x)   such  that  u(o,x)   Is  continuous  and  we 
usually  require  that   3u/^x   also  be  continuous.   We  choose 
the  Initial  values  of  p  and   p   such  that  the  solution  Is 
a  forward  facing  simple  wave  If   R~  -   0,      that  Is 


p(0,x)  =  p 


o 


1  +21zl  iL 
2  ^o 


7-1 


p(0,x)  =  -^  p(0,x)^  . 

We  fix  the  boundary  conditions  by  requiring  that  u(t  ,0)  =  u(0,0) 
and  u(t  ,c)  =  u(0,c)   for  all   n.   We  must  choose   c   large 
enough  so  that   c   remains  downstream  from  the  compression  or 
rarefaction  wave.   These  boundary  conditions  are  somewhat 
dubious  If   R~  7^  0  since  the  Navler-Stokes  equations  are 
not  hyperbolic,  but  the  results  appear  reasonable.   If  R~  =0 
and  u(0,0)  =  1  (u(0,0)  =  -  l),   then  we  have  a  compression 
(rarefaction)  wave.   In  the  former  case  a  shock  will  ultimately 
appear  in  the  flow.   In  the  latter  case  we  can  compute  the 
exact  solution   (R~   =  0,   thus  we  have  a  simple  wave)  by 
solving  a  simple  relation  involving  j^(x)  [4].   Thus  we  have 
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a  check  on  the  accuracy  of  our  finite  difference  methods . 
The  results  of  the  computations  are  described  below.   The 
computations  were  done  on  the  IBM  709O  at  the  New  York 
University  Computing  Center. 

4 .    An  Analysis  of  Stability, 

By  using  the  method  of  von  Neumann  we  shall  derive  a 

stability  criterion  for  the  Iteration  method  for  the  case 

where   R   =0.   To  use  this  method  we  must  assume  that 

f (v)  =  Av  where   A  Is  a  constant  matrix  (A  has  real  distinct 

eigenvalues  In  our  case)  [5].   Then  the  Iteration  method  Is 

given  by 

o    n 
w   =  V 

n+1    p 
V    =  w^ 


w 


-S+1  =  v'^  -  -^  (AwS  +  Ava)  =  v"^  +  Nv'^  +  Nw^   where   Nv  denotes 


■j^     ViiW^  i-  AV^ 


the  operator   -  -jj-  Av^  .   We  see  that 


V 


^^+1  =  wP  =  (I  +  2N  +  2n2  +...+  2NP)v'^ 


To  use  the  method  of  von  Neumann  we  assume  that   v.  =  K  exp(jwx^) 
where   j  =  \/^  .   Then  we  must  evaluate  the  eigenvalues  of  the 
amplification  matrix  M  where   v"^"^"^  =  Mv"^.   If  m-  =  -  (A/2)slnwAx, 
then  Nv"^  =  j^iAv^ .   Therefore   M  =  I  +  2[ijA  +•••+  2j'^\i^Pp .      If  we 
let  m  and  a  denote  the  eigenvalues  of  M  and  A,   and  let 
b  =  ixa,   then 

m=-l+2(l+jb+...+  jPbP), 
"1  -  ,-P+lbP+^ 


or     m  =  -  1  +  2 


1  -  Jb 
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Therefore  we  have 


\mf   =   1  +  (-D^^b^^-^+l) 


1  +  (-dV^'^+i^' 


and 


|m|2  =  1  +  (-D'^+iiibS^^+i) 


1  +  b^ 


1  +  (-D^+H^n 


If   p  =  2n  +  1, 


1  +  b' 


If   p  =  2n 


Prom  these  equations  we  see  that   |m|  >  1   if   p  =  1,  2,  5, 

6,    9,   10,  ,      or  If   jbl  >  1.   If  p  =  3,  ^,   7,   8,  

and   |b|  ^1,      then   |m|  ^  1   which  shows  that  the  method  Is 
stable  In  this  case  and  unstable  In  the  preceding  case.   Note 
that  the  Courant,  Prledrlchs,  Lewy  condition  Is   1^1  i  p-  • 

Computations  were  performed  using  the  Iteration  method 
as  described  above  In  section  5-   In  one  run  we  had 
0ix)    =    (x-d)  (2d  +  3  -  2x)  -  1   and  therefore  the  solution 
was  a  rarefaction  wave   (R   =  0  In  all  the  tests  of  stability 
for  the  Iteration  method).   With  b  =  0.95,  Ax  =  0.025,   and 
using   p  =  2   Instability  was  evident  after  39  time  steps, 
at  which  time  the  pressure  took  on  negative  values.   We  also  used 


(4.1) 


0(x)  =  (x-d)  (2x  -  2d  -  3)  +  1 

+  (l/5)(x-d)(d  +  1  -  x)sln([(l/4)Ax]7r(x-d)  ) 


([a]  =  Integer  part  of   a)   which  produced  a  compression  wave 
with  a  superimposed  oscillation.   With  b  =  0.95,  Ax  =  0.025 
and  p  =  2   Instability  occurred  at  10  time  steps  and  with 
p  =  5  at  l4  time  steps.   With  p  =  3?  4  or  7 
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both  problems  were  run  to  100  time-steps  with  no  sign  of 
Instability.   In  the  case  of  the  compression  wave  a  shock 
eventually  formed  in  the  flow.   This  did  not  cause  Instability, 

One  run  was  made  where  the  solution  was  a  rarefaction 
wave  with  b  =  0.45,  Ax  =  O.O5,   and  p  =  2.   No  instability 
was  evident  after  200  time  steps «   This  can  probably  be 
explained  by  inspection  of  the  eigenvalues   m  of  the  ampli- 
fication matrix  M  defined  above.   If  p  =  2,   then 
|m|  =  1  +  4b  .   Since   log(l  +  4-0.95)  =  O.63   and 
logd  +  4-0.45)  -   0.064,   we  see  that  approximately  ten  times 
as  many  cycles  will  be  required  to  produce  instability  when 
b  =  0.45.   Note  that  the  amplification  after  N   time  steps 
is  governed  by  m  ,   i.e.   (1  +  4b  )  . 

When  p  =  2   the  iteration  method  is  very  similar  to 
the  Lax-Wendroff  method.   For  p  =  2   and  R~   =0   the 
equations  for  the  iteration  method  are  (assume  that  f(v)  =  Av, 
where   A  is  constant) 

2 
/ 1,  ^^    n+1    n   A./n      n\,A.2/n     „n,n\ 
(4.2)    V.    =  V.  -  2  A(v.^^  -  v._^)  +  ^  A  (v.^2  -  2Vi  +  v.,^) 


and  the  Lax-Wendroff  equations  are 

2 
f  ],    -,\  n+1    n   A„/n      n\,A.2/-n     ^n,n\ 

(4.5)    V.    =  V.  -  2  A(v.^^  -  v._^)  +  2  A  (v.^^  -  2v.  +  v._^) 

The  Lax-Wendroff  equations  are  derived  from  equations  (2.1) 

by  replacing  derivatives  by  difference  quotients.   If  we 

n    _   n    n 

,   ^  V   ,     ^i+2      1    ^1-2    .   ,   ,   „/-,/,  x2   n 
replaced  — ^       by   ?5 instead  of  (1/Axj   v^-, 

bx^  4(Ax)2  '''' 
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then  the  Lax-Wendroff  method  would  yield  equations  (4.2) 
instead  of  (4.3)=   In  this  case  the  Lax-Wendroff  method  would 
produce  an  unstable  difference  scheme. 

We  are  unable  to  compute  a  stability  criterion  for  the 
Lax-Wendroff  method  as  applied  to  the  Navier-Stokes  equations. 
We  can  obtain  an  empirical  stability  criterion  by  computing 
the  solution  for  a  compression  wave,  for  various  values  of 
R,  Ax.     and  At.   We  compute  the  solution  for  100  time-steps 
and  observe  the  value  of  At   for  which  the  computation  becomes 
unstable  at  fixed  values  of  Ax  and  R.   In  these  computations 
we  use  the  Lax-Wendroff  method  on  the  equation   (2.2).   The 
initial  conditions  are  determined  by  0{x)      as  given  in  equation 
(4.1).   We  assume  the  scheme  is  stable  for  a  given  At   if 
there  is  no  obvious  instability  after  100  time  steps.   The 
results  of  the  computations  are  shown  In  table  1.   The  largest 
value  of  At   for  which  the  calculation  is  stable  is  given 
along  with  the  minimum  value  of  At  for  which  the  calculation 
is  unstable.   The  simple  explicit  scheme  for  the  heat  equation 
^u/bt   =  o^h   u/bx^      has  the  stability  criterion   At  =  ( l/2o)  (Z!oc  )^ . 
The  derivative  of  second  order  in  the  Navier-Stokes  equations 
(2.2)  has  the  coefficient   4/5Rp .   If  the  Navier-Stokes  equations 
behave  as  a  parabolic  system  we  might  therefore  expect  the  sta- 
blllty  criterion  At  =  (3Rp/8)(Ax)  .   If  the  behavior  is  hyper- 
bolic we  expect  the  criterion  At  =  (i/(u+c))Ax  where  the 
maximum  value  of  u  +  c   at   t  =  0   is  used.   From  table  1  we 

see  that  the  behavior  is  parabolic,  that  is   At-v(Ax)  ,   for 

4 
R  =  1   and  R  =  10.   It  is  hyperbolic,   At~Ax   for  R  =  10  .   . 
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As  far  as  these  results  are  concerned  we  see  that  the 
following  empirical  stability  criterion  will  suffice 


At  -  min 


[(Ax/(u+c)),  (l/2)(5Rp/8)(Ax)^ 


We  are  forced  to  halve  the  usual  parabolic  condition.   In 
table  1  the  range  of  At   in  which  instability  occurs  for 
various   Ax   and  R   is  shown  along  with  the  values  of  At 
given  by  the  various  stability  criteria. 

The  iteration  method  was  also  used  to  compute  solutions 
to  the  Navler-Stokes  equations.   The  values  obtained  for  the 
pressure  using  this  method  agreed  with  the  values  obtained 
using  the  Lax-Wendroff  method  (applied  to  equations  (2.2)), 
to  within  0.3°/o  for  Ax  =  O.O5,  R  =  10  at  40  time  steps. 
This  agreement  seems  to  be  consistent  with  the  accuracy  of 
these  methods  as  discussed  in  the  following  section. 

In  figure  1  the  values  of  velocity  obtained  using  the 
Lax-Wendroff  method  at   R~   =0,1   and  R~   =  0  are  plotted. 
Both  values  are  at  time   t  =  O.56,   All  initial  conditions 
are  the  same  for  both  calculations. 

5.    The  Accuracy  of  These  Finite  Difference  Methods. 

The  truncation  error  of  both  the  iteration  method  and 
the  Lax-Wendroff  method  is   0(Ax-^).   We  checked  this  only  for 
R"-^  =  Oj   it  is  undoubtedly  true  for  R""*"  7^  0  as  well.   If 
we  take   R~  =0  and  choose   0(x)   such  that  the  solution 
is  a  rarefaction  wave,  then  we  can  compute  the  exact  solution 
by  solving  an  equation  involving  0(x).   We  compare  this 
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solution  with  that  obtained  by  the  finite  difference  schemes. 
This  Is  done  In  tables   A,  B  and   C   below.   These  tables 
give  the  maximum  percentage  error  In  the  pressure.   The  error 
In  the  Lax-Wendroff  method  Is  not  greatly  different  from  that 
In  the  iteration  method. 

The  error  is  certainly  not   0(Ax-^^).   The  error  might  be 
considerably  reduced  if  the  exact  solution  were  analytic,  at 
least  smoother  than  the  ones  used  here.   This  is  determined 
by  the  choice  of  u(0,x).   In  the  case  of  tables   A  and  B 
the  second  derivative  of  the  exact  solution  is  not  continuous 

since  h0   /Sx  y^  0     at  x  =  d  and  x  =  1  +  d  when   t  =  0. 

1  p  p 

That  is   use   but  u  ^  c  .   In  table  C   u  e  c   but 

u  /  c^ .   These  difference  methods  have  been  used  on  a  problem 

for  which  the  solution  is  analytic  [?]•   In  this  case  the 

error  was   O(Ax^). 

Table  A.   The  maximum  percentage  error  in  the  values  of 

pressure  obtained  by  the  Iteration  method  when  u(t,x)  e  c  . 


Ax 

0.1 

0.05 

0.01 


Tlme-s 

teps 

20 

40 

100 

2.1 

1.9 

1.8 

0.66 

0.82 

0.90 

0.03^ 

0.056 

0.099 

Table  B.   The  maximum  percentage  error  in  the  values  of 
pressure  obtained  by  the  Lax-Wendroff  method  when  u(t,x)  e  c 
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Ax 

0.1 

0.05 

0.01 


Time-steps 
20       40 


100 


1.54 

1.4 

1.35 

0.48 

0.65 

0.76 

0.024 

0.041 

0.071 

Time-s 

taps 

20 

40 

100 

1.2 

1.33 

1.31 

0.36 

0.52 

0.65 

0.006 

0.015 

0.04 

Table  C.   The  maximum  percentage  error  In  the  values  of 
pressure  obtained  by  the  Lax-Wendroff  method  when  v(t,x)  s  c 

Ax 

0.1 

0.05 

0.01 

The  iteration  method,  using  three  iterations  per  time- 
step,  requires  240  seconds  to  compute  100  time-steps  when 
300  mesh  points  are  used  (an  IBM  7090  is  used  for  these  calcu- 
lations )  .   This  is   2.7  x  10  -^^   seconds  per  iteration  per  mesh 
point.   The  Lax-Wendroff  method  also  requires  240  seconds  to 
compute  100  time-steps  when  3OO  mesh  points  are  used.   Note 
that  both  these  calculations  are  for  the  Navier-Stokes  equations 
which  require  more  time  than  the  equations  of  inviscid  flow. 

6 .    Computation  for  Flows  Which  Contain  a  Shock. 

If  the  flow  contains  a  shock  these  finite  difference  methods 
can  still  be  used  to  compute  an  approximate  solution   [l,2]. 
The  different  methods  yield  results  which  differ  considerably. 
The  Lax-Wendroff  method  yields  much  better  results  when  the 
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equations  of  motion  are  written  in  conservation  form  as 

Professor  Lax  has  Indicated  [l,2].   The  iteration  method 

fails  even  though  the  equations  are  in  conservation  form. 

These  results  are  obtained  by  computing  the  solution 

for  a  compression  wave   (R   =  O) .   Enough  time  steps  are 

used  so  that  the  shock  created  by  the  compression  wave 

reaches  an  apparent  steady  state.   All  of  the  methods  seem 

to  be  stable  provided  the  shock  strength  is  not  too  great. 

In  figure  2  the  values  of  u   computed  by  the  iteration  method 

are  compared  with  those  computed  by  the  La_x-Wendroff  method 

(equations  in  conservation  form)  for  a  shock  of  strength 

((P-,  -  P  )/p  )  =  0.40.   Clearly  the  iteration  method  is  not 
1    o  '  "^o 

acceptable  even  for  this  weak  shock. 

In  figure  2  the  values  of  u   obtained  by  the  Lax-Wendroff 
method  (conservation  form)  are  shown.   In  table  2  the  state 
behind  the  shock  is  given  as  calculated  by  the  Lax-Wendroff 
method,  both  for  the  equations  In  conservation  form  and  in 
the  form  of  (2.2).   Both  the  minimum  values  and  the  values 
at  the  peak  immediately  behind  the  shock  are  given  (see  figure  2 

For  weak  shocks  we  can  obtain  an  approximate  value  for  the 
state  immediately  behind  the  shock  by  assuming  that  the  flow 
behind  the  shock  is  Isentroplc.   Then  we  may  assume  that  the 
characteristics  behind  the  shock  are  straight  lines.   Thus 
we  can  determine  the  value  of  u  +  c   behind  the  shock.   The 
state  behind  the  shock  then  follows  from  the  Hugoniot  relations. 
This  approximate  solution  is  given  in  the  tables.   The  values 
for  the  constant  state  behind  the  compression  wave  are  also 
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given.   These  are  the  values  used  to  determine  the  left-hand 
boundary  condition. 

If  we  know  the  pressure  jump  across  a  shock  the  velocity 
and  density  can  be  obtained  from  the  Hugonlot  relations.   Thus 
we  can  use  the  pressure  calculated  by  the  difference  scheme  to 
obtain  a  check  on  the  values  of   u  and   p.   The  values  of  u 
and   p   computed  from  the  Hugonlot  relations  are  shown  In  the 
table  . 

We  see  from  table  2  that  the  values  computed  with  the 
equations  in  conservation  form  more  nearly  satisfy  the  Hugonlot 
relations.   Also,  the  method  based  on  equations  (2.2)  Is  unstable 
at  a  lower  shock  strength  than  the  method  based  on  the  equations 
In  conservation  form.   The  latter  method  Is  unstable  when  the 
shock  strength  equals  1J>}    and  perhaps  for  shocks  of  lower  strength 
although  It  seems  to  be  stable  when  the  shock  strength  equals  6. 
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Ax 

10^  At 
(stable) 

10^  At 
(unstable ) 

10^(Ax/u+c) 

10^(1/2) (3Rp/8)Ax2 

0.05 

0 .  056 

0.061 

2.1 

0.04? 

0.025 

0.014 

0.015 

1.0 

0.012 

0.01 

0.0022 

0.0024 

0.42 

0.0019 

0.05 

0.56 

0.61 

2.1 

0.47 

0.025 

0.l4 

0.15 

1.0 

0.12 

0.01 

0.022 

0.024 

0.42 

0.019 

0.05 

2.0 

2.1 

2.1 

4.7 

0.025 

1.0 

1.1 

1.0 

1.2 

0.01 

0.22 

0.24 

0.42 

0.19 

0.05 

2.1 

2.3 

2.1 

470. 

0.025 

1.0 

1.1 

1.0 

120. 

0.01 

0.40 

0.43 

Table   1 

0.42 

19. 
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Const .state 

Approx .solu . 

Calcul 

ated 

From  Hug 

;oniot 

behind  the 

(isentropic ) 

solut 

ion 

relations 

shock 

min . 

max . 

min. 

max. 

p    1.49 

1.48 

1.48 

1.51 

1.48 

1.52 

Method  based 

P    7-20 

7.17 

7.15 

7.48 

-- 

-- 

on  equations 
in  conser- 

u   1.00 

0.995 

0.99 

1.07 

.99 

1.08 

vation  form 

p    2.l6 

2.04 

2.06 

2.11 

2.06 

2.12 

P    3.02 

2.92 

2.99 

3.10 

-- 

-- 

II        II 

u    1.00 

0.984 

1.00 

1.04 

1.00 

1.04 

p     4.21 

3.13 

3.20 

3.27 

3.31 

3.33 

P    1.93 

1.58 

1.80 

1.83 

-- 

-- 

II       (1 

u    1.00 

0.955 

1.04 

1.05 

1.04 

1.05 

P    1.49 

1.48 

1.47 

1.50 

1.46 

1.51 

Method  based 

P    7.20 

7.17 

7.04 

7.37 

-- 

-- 

on  equations 
(2.2) 

u    1.00 

0.995 

1.03 

1.13 

0.96 

1.05 

p     2.l6 

2.04 

2.02 

2.13 

1.96 

2.05 

P     3-02 

2.92 

2.74 

2.92 

-- 

__ 

M              11 

u     1.00 

0.984 

1.09 

1.15 

0.92 

0.97 

p     4.21 

3.13 

P    1.93 

1.58 

unstabl 

e 

II               II 

u    1.00 

0 . 9  55 

Table  2.   The  state  immediately  behind  the  shock  as  calculated  by 
the  Lax-Wendroff  method. 
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